Casimir interaction between topological insulators with finite surface band gap 
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Casimir interaction between topological insulators with opposite topological magnetoelectric po- 
larizabilities and finite surface band gaps has been investigated. For large surface band gap 
limit(m — > oo), we can obtain results given in [Phys. Rev. Lett. 106, 020403 (2011)]. For 
small surface band gap limit (m — > 0), Casimir interaction between topological insulators is attrac- 
tive and analogy to ideal mental in short separation limit. Generally, there is a critical value m c 
and when the surface band gap is greater than the critical value, the Casimir force is repulsive in 
an intermediate separation region. We estimate the critical surface band gap m c ~ l/(2a), where a 
is a critical separation where Casimir force vanishes. 

PACS numbers: 12.20. Ds, 41.20.-q, 73.20.-r 



1. INTRODUCTION 

Time reversal invariant topological insulator(TI) 1-3 is 
a new quantum state of matter which has a full insulating 
gap in the bulk, but has gapless surface states protected 
topologically. This material has been extensively stud- 
ied experimentally 4 ~ ; ' and theoretically 10-15 . Two dimen- 
sional TI has been observed in HgTe quantum well ,15 , 
Sbi-^Tea; is the first material has been reported to be 
3-dimcnsional TI, and Bi2Se3, Bi2Te3, Sb2Te3 have been 
predicted 11 ' to be TI with single Dirac cone on the sur- 
face. Novel properties of TI have been predicted, for 
instance, effective monopole ' and topological magneto- 
electric effect 13 , superconductor proximity effect induced 
Majorana fermion states 1 * etc. 

Recently, an interesting property of TI, tunable re- 
pulsive Casimir interaction between TIs with opposite 
topological magnetoelectric polarizability has been 
proposed 11 ', and the robustness of this repulsion in small 
separation limit against finite temperature and uniaxial 
anisotropy has also been analyzed . Repulsive Casimir 
interaction has been discussed in a few proposals, with 
special geometry 21 or chiral metamaterials 22 , or filling 
high-refractive liquid between dielectrics" . The repul- 
sion between TIs is analogy to metamaterials, however, 
time reversal invariant TI is protected by gapless surface 
states. In order to observe the repulsive Casimir interac- 
tion, one need cover the TI surfaces with magnetic coat- 
ing to open the band gap. The effect of finite surface 
band gap on this repulsive force is considerable. 

In this paper, we analyze the influence of finite sur- 
face band gap on Casimir force between TIs with oppo- 
site topological magnetoelectric polarizability 0, we show 
that there is a minimal surface band gap m c and when 
surface band gap m < m c , repulsive Casimir force will 
disappear. We also estimate this critical surface band 



gap numerically. 

Let us formulate the model. When time reversal sym- 
metry is protected in the bulk, the topological nontrivial 
term a/(47r 2 ) / cPxdtBE ■ B can be reexpressed as spin- 
momentum locked fermions on the interface of TI and 
normal insulator, in this paper we consider only one kind 
of fermion corresponding to = ir or — 7T, generaliza- 
tion to mult i- fermions is straightforward. Action of Dirac 
fermion on TI surface is 
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d 3 x ip [ij a (d a + ieA a ) - m] ip, 



(1) 
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where a = 0,x,y, 7 — a z , j x 
a x,y,z are p au ij ma trices of spin, vf is the Fermi velocity 
of surface fermion, which has a magnitude of 10 -3 speed 
of light (we set h = c = 1 in this paper) and takes different 
values for different materials 5,1 ', ie, vf = 1.3 x 10~ 3 for 
Bi2Tc3, and 1.7 x 10 -3 for Bi2Sc3. Parameter m is sur- 
face band gap opened by magnetic coating on TI and we 
assume chemical potential has been tuned into the sur- 
face band gap. A a present the first three components of 
vector-potential, while electromagnetic field is described 
by Maxwell action: 



Sem = ~f d A x{sE 2 - -B 2 ), 
8ttJ fi 



(2) 
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where E and B are electric and magnetic fields, e and fi 
are permittivity and permeability of TI in the bulk and 
equal to 1 in the vacuum. 

This paper is organized as follows: In Sec. 2, we eval- 
uate an effective action for electromagnetic field on TI 
surface by quantum field theory approach and give the 
Maxwell equations of electromagnetic field with proper 
boundary conditions. In Sec. 3, we analyze the Casimir 
interaction between TIs via Lifshitz theory. We discuss 
the results in Sec. 4, and give a conclusion in Sec. 5. 



2 



2. EFFECTIVE LAGRANGIAN ON TI 
SURFACE AND MAXWELL EQUATIONS 



In order to calculate the Casimir interaction caused 
by quantum fluctuation of electromagnetic field between 
TIs, one need to integrate the contribution from sur- 
face fermion. An effective action for external electro- 
magnetic field in (2+l)-dimension can be found by stan- 
dard quantum field theory approach , S e ft(A) = 
— i lndet[«7 a (<9 a + ieA a ) — to]. We introduce a Feynman 
parameter, integrate out the fermion field up to one-loop 
correction and get the effective action in the following 
form: 



Seff(A) 



d 3 x 



8tt 



e abc A a d b A c 



3(A) 
47t|to| 



F 0j F°i + v%F xv F xy 



F 1 xyi 



(3) 
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FIG. 1: Schematic illustration of Casimir interaction between 
TIs with opposite topological magnetoelectric polarizability 9. 
We assume the thickness of magnetic coating is much smaller 
than the separation between TIs. 



tions with surface corrections: 



with dimensionless parameters </> and $ which take the 
forms: 



f 1 

0(A) = sign(m)a / dx 

Jo \/l - x(l — x)\ 

$(A) = a [ dx . ^~ x "> x 
Jo ^/T=xJi 



y/l -x(l -X)\' 



(4) 
(5) 
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4tt 
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V D 

[d t D - (V x H)] 



5(z - Zi) 



47T 2tt to 



(8 t E - v 2 F V x B) 



V • B = 0, 

d t B + (Vx£) = 0, 



,(9) 

(10) 
(11) 



where sign(m) gives the sign of surface band gap, 
which corresponding to different signs of topological 
magnetoelectric polarizability. a = 1/137 is the fine 
structure constant of electromagnetic interaction, A = 
[kg — v F (k x + ky)] /to 2 , and fco, k x , k y are frequency 
and momentum of electromagnetic fields on TI surface. 
A detailed derivation and a short discussion on this ef- 
fective action(3) have been given in the appendix. Wc 
also note that in both limit, to 2 — > and to 2 — > oo, (ft 
and $ are convergent. For the sake of Eq.(20), we derive 
expressions of <p arid $ in imaginary time formalism: 



~i, \ . , ,2a / -v/7 i 

0(7) = sign(TO)-^ arctan [ ] , 



— + ( — Pttt ) arctan ( 

27 UV7 7 3/ V V 2 



(0) 
(7) 



where 7 = {k^ + v F {k 2 + k 2 ))/m 2 . For the large surface 

band gap limit (|to| — > 00), 4>{l) ~~ ^ sign(m)o:, the term 
proportional to (f>(\) in Eq.(3) is topological and the term 
proportional to $(A) in Eq.(3) is vanishing. For the small 
gap limit (|to| — > 0), 0(7) ~> and $(7) 1/6. 

Add the surface term Eq.(3) to standard action of elec- 
tromagnetic fields Eq.(2), one can get the Maxwell equa- 



where D = eE and H ~ B/fi are electric displacement 
field and magnetizing field, Ej = tj^Ek {j,k = x,y), 
(i = 1,2), z\ = and z-i = a arc positions of Tl-surfaccs 
(as shown in Fig.[l]), 0i and 02 are corresponding values 
of 0. Without loss of generality, we assume the absolute 
values of surface band gaps on two TIs are equal, differ- 
ent signs of surface band gaps stand for different signs 
of the topological term a9E ■ B^Att 2 ) in Lagrangian of 
electromagnetic fields in TIs. We also note that in large 
band gap limit (|to| — > 00), these Maxwell equations are 
equal to those given in Rcfs. 1 1 ' Ti by redefine the electric 
displacement and magnetizing field as D = eE + a^B, 
H = j^B — a^E. From above Maxwell equations, we 
get the following discontinuous boundary conditions: 



2<I> 



-<t>iB z + — (d x E x + d y Eyll2) 



D z (z+) - D z (zr) 

H x (z+) - H x (zr) = <j>iE x -^(d t Ey + v F d x B z ()13) 
Hy(zt) - Hy(z-) 



2$ 



24> 



E y + -— (d t E x - v 2 F d y Bj)\A) 



where z\ means z^ 
on the interfaces. 



± 0. And E x , E v , B z are continuous 
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CASIMIR INTERACTION 



Now we analyze the Fresncl coefficients of reflection 
light on the Tl-vacuum interface. Incident TE-modc 
from vacuum with wave- vector (k x ,k y ,k z ) will induce re- 
flected TE and TM-mode, we assume the reflection co- 
efficients are r ee and r em respectively, then the electro- 
magnetic waves in the vacuum read: 



E = (1 + r ee )k (-k y e x + k x e y ) 
B = {-k z k + k 2 e z )+r ee (k z k~\ 

~l~^'em^o( kyG x -\- h x Gy)^ 



i{ k z h k 62), 



k 2 e z ) 



(15) 



and the refracted light with TE, TM-mode in TI take the 
forms: 

E = t ee k (-k y e x + k x e y ) + ct em (p z k - k 2 e z ), 
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t ee (-Pzk + k 2 e z 



-k Q (-k y e x + k x e y )(16) 



where t ee and t em are refraction coefficients of TE and 
TM-mode, c is the relative velocity of light in TI bulk, 
k = k x e x + k y e y , k 2 = k 2 + k 2 and p z is z-component 
of wave vector in TI. For the injected TM-mode, one 
can write the analogy equations with reflection coeffi- 
cients r me , r mm and refraction coefficients t me , t mm . Af- 
ter some tedious derivation, we obtain the Fresnel coeffi- 
cients matrix 1Z in imaginary time formalism: 



with 



K 



= -1 



' me ' mm 



2 
D 



2 
D 
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1 + e- 



—2i 

fi k z 



2$^ 
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2A$— 

k z 



(17) 
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FIG. 2: The ratio E^/E^ as a function of dimensionless 
separation auj for different oscillator strength g' = ^/gj/coj 
in the closed surface band gap limit, m = 0. Where 
Eq (Eq ) is the Casimir energy with(without) surface cor- 
rection. Here fermion velocity vf = 1.0 x 10 -3 . 

where A is the surface area of TIs, TC- 12 ' are Fresnel 
coefficients on the surfaces, k$ - 



fen ~\~ k'fl . 



In order to 



calculate the Casimir energy density numerically, we also 
need a form of frequency-dependent dielectric permittiv- 
ity e (we assume the permeability fj, = 1), this can be 
modeled by 28 ' 29 : 



e(iko) = 1 + 



K 



9j 



(21) 



we consider only one oscillator (K = 1) with oscillator 
strength gj, oscillator frequency ujj and damping param- 
eter 77. 7,/ <C loj and we omit the contribution from 
damping parameter here. 



where the denominator 
D = [ 



1 



£—11 

2$^ 

m 



2 7 $ — 

Kr. 



1 Vr 

fx k z 



1 _4 7 $ 2 ) .(19) 



For the large surface band gap limit, we can obtain the 
same Kerr rotation and Faraday rotation angle as given 
in Ref. 13 ' 27 . 

In imaginary time formalism, Casimir energy density 
between two parallel dielectric semispaces can be ex- 
pressed in a closed form of dielectric permittivity: 



E c (a) 
A 



dko 
"27 



d 2 k 



(2tt)= 



log dct 



n (l) n (2) e -2k 3 a 



(20) 



RESULTS AND DISCUSSION 



In this paper, we analyze 
TIs with finite surface band 
band gap limit (m — > 00), 
given in [Phys. Rev. Lett, 
equations (17)-(21). Second 
limit (m —> 0), the off-dia^ 
ficients matrices will vanish 
rewritten in imaginary time 



Casimir interaction between 
gap. First, for large surface 
we can obtain same results 
106, 020403 (2011)] 1!l from 
, for small surface band gap 
;onal terms in Fresnel coef- 
and Casimir energy can be 
formalism as: 



E C (a) 



dko 
log (l 



d 2 k\\ 

-2k 



log(l 



-2ka (1) 

TE TE ) 



Jl) r (2) 
'TM'TM 



(2)\ 
~E J 

(22) 
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with 



r TE = -1 







e + 




2 ff 


\J cos 2 



Up, sin 



D p Sill 



7TQ: £3 

4 fc 3 



where fa 



fc 



fen 



u 3 



P3 



til sin 2 i 



= , (23) 



= (24) 



fc 2 , and = 



cos 1 (k /k 3 ). 

The Casimir energy between dielectric materials with- 
out special boundary conditions, a — > in Eq.(23) and 
Eq.(24), has been studied 28-30 . 

Considering correction from surface interaction, for 
large separation limit, we obtain the correction up to 
first order of fine structure constant: 



lot 



not 
Ad? 



e(0) - 1 



(e(0) 
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FIG. 3: Boundary of repulsive and attractive Casimir inter- 
action in the plane of dimensionless separation d — atuj and 
product \m\a for different oscillator strengths g' — ^fgj/uij. 
When the parameters d and \m\a have been taken in the left- 
up region over these lines, the Casimir interaction is attrac- 
tive, when d and \m\a have been taken in the lower-right re- 
gion of these lines, the Casimir interaction is repulsive. (The 
relative Fermi velocity up has been taken to be f.O x 10 -3 ). 



where Eq = Alu 3 j / (2ir) 2 which is set as the unit of Casimir 
energy, d = auij is the dimensionless separation. 

For small separation limit, in order to make the physics 
more clear, we also formally expand Eq.(22) in powers of 
a, up to first order correction, the Casimir energy takes 
the following form(here we assume the relative oscillator 
4 



strength gj/ojj < 1): 



E { c\a) 
Eq 



g^n^ f°° 2 

^64^ dy y e 

~8(t) r 6(-t) 
— ^arctanvi H ==-axctanh\ 



,(26) 



where t = —1 + vpy 2 /Ad 2 and 6 if) is the Heaviside unit 
step function. Casimir energy is dominated by surface 
Dirac fermion and turns into the ideal conductor case 
which is proportional to 1/a 3 . This conclusion is also 
confirmed numerically in Fig. [2]. 

Finally, for general surface band gap, we have two di- 
mensionless parameters: m/uj and d = uja (there are 
two other parameters in our model, the Fermi velocity 
of surface fermion, Vf, and optical oscillator strength in 
TIs, gj/uij, which both have quantitatively influence on 
Casimir energy). For the large separation limit (a 3> 
max(l/cjj, l/|m|)), we expand the integral in Eq.(20) in 
power of fine structure constant 51 a and consider the cor- 
rection up to a. In this case, the dielectric permittivity 
s{iko) can be approximated by long wave length limit 
value e(0), and the Casimir energy correction from in- 
teraction between surface fermions and electromagnetic 
field reads: 



E { c\a) 



En 



\m\a 
ujd? 



dx 



r 2 (e(0)-r) x 2 (l-r)' 
(e(0) +r) 3 + (1 + r) 3 
(27)" 

where r = \Jl + (e(0) — l)a; 2 , and vf <C 1. 

For the small separation limit (d — > 0), in order to 
make the physics more clear, we also formally expand 
the Casimir energy in power series of a. In this case, 
the Casimir energy is dominated by surface terms, the 
term which contains $ 2 and is proportional to l/m 2 a 5 
is important. However this dominant term will be sup- 
pressed if \m\a — > oo, the topological term which contains 
sign^i^)^ 2 and is proportional to 1/a 3 will provide a 
large repulsive potential between TIs when sign($i#2) = 
— 1. So the surface terms in Casimir energy will dominate 
and \m\a is a good parameter to estimate the Casimir 
force: when \m\a S> 1, the Casimir force will be repulsive 
and when \m\a <C 1, the Casimir force will be attractive. 

In Fig. [3], we give the boundary of repulsive and at- 
tractive Casimir interaction, as a function of dimen- 
sionless separation d = auij and product \m\a. We 
find that there is a critical value (|m|a) c ~ 1/2, when 
\m\a < (\m\a) c , the Casimir interaction is attractive for 
any separation length. The independence of (|m|a) c on 
oscillator strength gj /lo 2 shows that Casimir interaction, 
in small separation limit, is dominated by surface terms. 
More intuitively, we calculated the Casimir energy as a 
function of dimensionless separation d = aujj for different 
surface band gaps, as shown in Fig. [4], for given param- 
eters, g,j/uj 2 = 0.45 2 and v F = 1.0 x 10~ 3 . We find 
the critical surface band gap, where the repulsive peak 
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FIG. 4: Casimir energy density _Ec (in units of Eq = 
ojj/(2ty) 2 ) as a function of the dimensionless distance d — aujj 
with different surface band gaps m/ujj. Here we take the 
dimensionless oscillator strength g.;/to^ = 0.45 2 and Fermi 
velocity v F = 1.0 x 10" 3 . 



FIG. 5: Boundary of repulsive and attractive Casimir inter- 
action in the plane of dimensionless separation d — aujj and 
product \m\a for different topological magnetoelectric polariz- 
ability. (The relative Fermi velocity up and oscillator strength 
gj/uij has been taken to bel.OxlCP 3 and 0.45 2 respectively). 



disappears, m c 300wj (the blue-square dotted line in 
Fig. [4]). 

We note that our calculations can be generalized to 
multi-value of topological magnetoelectric polarizability 
9 = (2n + l)7r (n is an integer) straightforward by in- 
troducing multi-fermion on TI surface, and the critical 
value (|m|a) c is independent on the absolute value of 9 (as 
shown in Fig. [5]), this is because in short separation limit, 
Casimir interaction is dominated by surface terms and 
each species fermion will contribute both repulsive and 
attractive Casimir interaction if sign(#i) = — sign(#2)- 

We can use this relationship to estimate the critical 
surface band gap for repulsive Casimir interaction. For 
TlBiSc2 suggested in Ref. 10 , the minimum of Casimir 
energy appears at a separation of a ~ 0.1/xm, and the 
corresponding surface band gap needs to be greater than 
leV, which reflects that the width of surface band gap 
opened by magnetic coating is non-ignorable and unac- 
cessible experimentally. 



5. CONCLUSION 

We studied the Casimir energy between TIs with oppo- 
site topological magnetoelectric polarizability and finite 
surface band gap via Lifshitz formula, we found that, 
in small separation limit, Casimir force is dominated 
by interaction between surface fermion and electromag- 
netic field in the vacuum, and a great surface band gap 
m > m c ~ l/(2a) is essential for repulsive Casimir inter- 
action. 
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Appendix: Effective action 

We give a detailed derivation of the effective action(3) 
in this appendix. The effective action from quantum field 
theory is: 

Seff(A) = -J — Aa (k)n ab (k)A b (k), (A.l) 

where n(fc) is the polarization tensor, which takes the 
form: 

iIP 6 (fc) = -e 2 j ^tr[H 7 a )G(fc+p)H 7 b )G(fc)], 

(A.2) 

and G(k) = i/ {^ a k a + m) is the propagator of fermion on 
Tl surface. From the standard calculation in quantum 
field theory, one can get the exact form of polarization 
tensor: 

n(fc) = rii(fc) + n 2 (fc) (a.3) 
nf (k) = t^ e abc ik , A4) 

4-7T 

<j)(Al / k x + k y —kok x —k()ky \ 

n 2 (fc) = I -k k x kl-v 2 F k 2 y v 2 F k x k v J 

\ k()ky Vpk x ky Vph x J 

(A.5) 

where 0(A) and $(A) has been given in Eq.(4) and Eq.(5), 
^1,2(^0) are the momentum(frcqucncy) of electromag- 
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nctic field. One can check that the polarization tensor 
satisfies Ward identity, J2 a k a^ ab (k) = £ b n a6 (fc)fc fc = 0. 
The Fourier transformation of Eq.(A.l) gives Eq.(3). 



We take W y (k) as an example to show more detailed 
calculations of polarization tensor. Taking the trace in 
Eq.(A.2), one can get 



J 



iW y {k) = 



d 3 p 2vp[—imko + v F (2k x k y + k x p y + k y p x )] 



(271)3 [(p 0+ fc )2+ m 2 



r 



v 2 F (p- 



(A.6) 



One can get the following form of iW v (k) by introducing 
a Feynman parameter x and redefining the integration 
variables l' a = p a + xk a , Iq = 1' , and I = vpl': 



i\\ xy {k) = -2e 2 



dx 



d 3 l imko + 2x(l — x)v F k x k y 



(2tt) 3 



(A.7) 

where A = m 2 — x(l — x)(fcg — v F k 2 ) = m 2 [l — Xx(l — x)]. 
Making Wick rotation Iq — > ily and integration over I, 
we find: 



iW v {k) 



dx 



+ Vpk x k y 



imko 

J 4 

47T 2tt m 



dx 



x(l — x) 



27r\/A 
(A. 



Comparing with the effective action of electromagnetic 
field in monolayer graphene system as shown in Ref/'~, 
we find that there is an additional topological term 
Eq.(A.4) together with the normal vacuum polarization 
Eq.(A.5), the first term is essential for TI because 
this parity-odd term reflects the fact that there are 
always odd species of surface fermions which are spin- 
momentum locked, the contribution from second term is 
analogy to Dirac fermion in monolayer graphene system 
and reflects the dynamical response of TI surface state 
to extra electromagnetic field. 
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